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Abstract 

The electronic band structures of orthorhombic (oP28) and monoclinic (mC28) MnSb2S*4 were 
investigated with ab initio calculations in the local spin density approximation (LSDA) to the den- 
sity functional theory (DFT). An analysis of the electronic properties and of the chemical bonding 
is provided using the augmented spherical wave (ASW) method considering nonmagnetic, ferro- 
magnetic, ferrimagnetic and antiferromagnetic model orderings. In agreement with experimental 
results both modifications of MnSb2S4 are predicted to be antiferromagnetic. While the experi- 
mental band gap is missed for the monoclinic polymorph, the calculated band gap for orthorhombic 
MnSb2S4 is close to the experimental one. 
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I. INTRODUCTION 



Magnetic and semiconducting manganese sulphides attracted attention of solid 
state chemists since the early days of X-ray crystallography and magnetic structure 
investigations 1 " 5 done on haurite (MnS 2 ) and alabandite (MnS). They exhibit high magnetic 
moments due to the coordination of Mn 2+ in MnS6 octahedra where it prefers a high spin 
state with five unpaired electrons. However, MnS2 has been discussed as a rare example of 
a high spin to low spin transition under high pressure 6,7 . 

In the past few years the chemistry of magnetic manganese materials was enriched by 
fascinating discoveries mainly on multinary manganese oxides 6-8 . Properties like the giant 
and colossal magnetoresistance (GMR, CMR) inspired new fields of research on magnetic 
semiconductors. Besides promising technological applications and experimental challenges 
there is an increasing demand and success of theoretical understanding of the underlying 
chemical bonding and electronic properties. The development and application of effective 
density functional (DFT) methods within the local spin density approximation (LSDA) still 
plays an increasing role herein 9-11 . 

Fascinating properties were also discovered on manganese chalcogenides. MnS and MnS 2 
show antiferromagnetic ordering while diluted magnetic semiconductors (DMS) based on 
MnS exhibit outstanding properties related to spintronic applications 12-14 . Multinary ma- 
terials like MnCr 2 S 4 provide additional potential with respect to anisotropic resistivity and 
magnetic properties 15,16 . 

Due to its reduced dimensionality MnSb 2 S 4 serves as a promising low dimensional mag- 
netic semiconducting material. Contrary to spinel type MnCr 2 S 4 with Mn 2+ in MnS 4 tetra- 
hedra one finds MnS 6 octahedra in MnSb 2 S 4 . Therein, it is related to MnS and MnS 2 as 
well as in the observation of phase transitions. Orthorhombic MnSb 2 S 4 is accessible by 
hydrothermal synthesis and was earlier shown 17 to be isotypic to FeSb 2 S 4 which is an an- 
tiferromagnetic material 18 . Recently a new monoclinic modification (mC28) of MnSb 2 S 4 
was synthesized by high temperature methods 19 . MnSb 2 S 4 (mC28) can be transformed 
reversibly into the orthorhombic modification (oP28) at high pressure 20 . By electrical con- 
ductivity and magnetic susceptibility measurements it was found that MnSb 2 S 4 (mC28) is 
a semiconducting antiferromagnet with = 26.5 K and an electronic band gap of 0.77 
eV 19,20 . Concerning the bonding situation one faces ID magnetic interactions, as well as 
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bonds with and within the [S'&Sy 3- ligand network that is related to Sb 2 S 3 24 . However, no 
theoretical investigations are reported yet. Considering MnS and MnS2 again as prominent 
examples, LSDA calculations 21-23 achieved good agreement with experimental results, i.e. 
the prediction of semiconducting and magnetic ground states with moments around 4.5 Hb 
for Mn 2+ . For a-MnS the antiferromagnetic ground state was correctly found 21,22 . LSDA 
total energy calculations on MnS2 supported the possibility of a low spin/high spin phase 
transition for a compressed cell 22 . 

To discuss the differences and relations of the bonding, spin states and magnetic or- 
dering in mC28 and oP28 MnSb 2 S4 first principles calculations are subsequently reported 
modelling non magnetic (NM), ferromagnetic (FM), ferrimagnetic (FIM for the monoclinic 
system) and antiferromagnetic (AFM) structures in order to identify the ground state con- 
figuration. The applied augmented spherical wave (ASW) method was successfully used 
in previous calculations on magnetic semiconducting manganites 11 . The crystal structures, 
computational details, and results of the calculations on non-magnetic, ferromagnetic, and 
antiferromagnetic configurations are presented, as well as electronic band structures, site 
projected densities of states and chemical bonding characteristics. 

II. CRYSTAL STRUCTURES OF ORTHORHOMBIC AND MONOCLINIC 
PHASES 

For the calculations presented herein, the crystal structures of both MnSb 2 S 4 modifica- 
tions as determined by single crystal X-ray diffraction were taken as the starting points 17,19,20 . 
The space groups and the relevant lattice parameters used in the calculation are given in the 
first part of Table 1. Both modifications are based on chains of edge-sharing MnS6 octahedra 
(Fig. 1). These chains of octahedra are linked by [SfrSy 3- units to form layers in the case 
of MnSb 2 S 4 (mC28) and a three dimensional (3D) network in the case of MnSb 2 S 4 (oP28). 

The Sb-S bonds determine both the structural anisotropics and the differences between 
the modifications. Sb atoms exhibit a 3+2+x (x — 1, 2) coordination with three Sb-S bonds 
of about 2.5 A and two Sb-S bonds between 2.9 and 3.1 A ("secondary bonds"). In addition, 
there are so-called non bonding distances 3.1 A< d(Sb-S) < 4 A. Distinguishing these three 
types of Sb-S interactions we find all short Sb-S bonds linking edge sharing MnSg-octahedra 
of one chain in MnSb 2 S 4 (mC28). Slightly longer bonds link the octahedra to form a layered 



3 



structure (Fig. 1). Between the layers (along the c axis) only so-called non-bonding Sb-S 
distances are found. In the case of MnSb 2 S 4 (oP28) one finds double chains of octahedra 
which are interlinked by short and secondary Sb-S bonds. These double chains form a kind of 
fishbone scheme and non-bonding Sb-S distances between them result in a 3D network. The 
density of the title compound increases from 4.24 g/cm 3 (mC28) 19 to 4.51 g/cm 3 (oP28) 
17 , showing that the orthorhombic modification is the high pressure form. The distances 
rf(Mn-S) vary from around 2.6 A (Table 1) in both modifications. Thus, they show a 
slightly broader range than in the pure manganese sulphides with octahedral coordination 
of manganese, i.e., d(Mn-S) = 2.61 A in a-MnS 2 and d(Mn-S) = 2.59 A in MnS 2 1 . There 
are two different Mn positions in MnSb2S4 (mC28), with a higher site symmetry than the 
single Mn position in MnSb 2 S4 (oP28). The distortions of the MnSg octahedra are due 
to the structural anisotropy imposed by the [S'frS's] 3- units. They result in tetragonally 
distorted MnS 6 octahedra with a coordination number of 2+4 in (mC28) MnSb 2 S 4 and a 
coordination number of 1+1+2+2 in MnSb 2 S 4 (oP28), respectively. Further details are 
provided in refs. 17 ' 19 . 

Considering the magnetic coupling of manganese in the two polymorphs of MnSb 2 S4 the 
structural anisotropy provided by the MnSg chains has to be kept in mind. Thus, only 
two contacts rf(Mn-Mn) 3.8 A are present in the title compound, and all other distances 
between Mn atoms are larger than 6 A. This situation is quite different from the isotropic 
Mn sublattices of, e.g., a-MnS (12 x d(Mn-Mn) « 3.7 A), and MnS 2 (12 x rf(Mn-Mn) « 4.3 
A). This allows to investigate FM models with equal Mn spin directions and AFM models 
with alternating Mn spins along the chains in MnSb 2 S4. 

III. COMPUTATIONAL FRAMEWORK 

The electronic properties have been self-consistently calculated in the framework of the 
density functional theory DFT 25 ' 26 using the ASW method as implemented by Williams et 
al. 27 and Eyert 28 . The effects of exchange and correlation were parameterized according 
to the local spin density approximation (LSDA) scheme of Vosko, Wilk and Nusair 29 . All 
valence electrons, including 4rf(Sb) ones, were treated as band states. In the minimum 
ASW basis set, we chose the outermost shells to represent the valence states and the matrix 
elements were constructed using partial waves up to l max . = 2 quantum number. The ASW 
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method uses the atomic sphere approximation (ASA) which assumes overlapping spheres 
centered on the atomic sites where the potential has a spherical symmetry. In order to 
represent the correct shape of the crystal potential in the large voids of the respective 
crystal structures, additional augmentation spheres were inserted 28 to avoid an otherwise 
too large overlap between the actual atomic spheres. 

The calculations implicit of zero entropy (T = K) were started assuming a non-magnetic 
configuration which is non spin polarized (NSP) meaning that spin degeneracy was enforced 
for all species (atoms and empty spheres). Note that this configuration does not translate a 
paramagnetic state which would actually require a supercell with different orientations of the 
spins over the crystal sites. In a second step spin polarized (SP) calculations were performed 
by initially allowing for differing spin occupations, i.e., majority (spin up f) an d minority 
(spin down |) spins for all atomic species. The occupancies were self-consistently changed 
until convergence of the total energy (AE< 10~ 6 Ry.) and of the charges (AQ< 1CT 6 ) 
between two subsequent iterations was reached. For that a sufficiently large number of k 
points was used with respect to selfconsistancy of the results. In view of the large cells 
especially when symmetry is broken by introducing the antiferromagnetic orderings, we 
used up to 12*12*12, i.e. 1728 points to produce respectively 216 and 468 k points in 
the irreducible wedges of the orthorhombic and monoclinic Brillouin zones. Calculations 
are implicit of collinear magnetic structures. However non-collinear magnetic structures 
can occur in manganese based compounds such as in the nitride Mn 4 N which was studied 
in the same calculational framework 30 . In fact such heavy calculations could be achieved 
with great accuracy in energy differences between the magnetic configurations provided one 
considers high symmetry structures such as that of cubic anti-perovskite Mn 4 N. When one 
magnetic/crystallographic sublattice of all species is accounted for, a ferromagnetic order 
(FM) is described. Two magnetic sublattices need to be accounted for to calculate the 
AFM configurations. This can be achieved by symmetry breaking of the system, half of the 
constituents being "spin up" and the other half being "spin down" . This approach accounts 
for the effect of low spin and high spin Mn 2+ and spin spin interactions in AFM and FM 
models for MnSb2S4 similar to the incommensurate magnetic structure of FeSb 2 S4 18 . Indeed, 
spin reorientation, spin disorder, and the competition between AFM and FM orientation are 
discussed to play an important role in magnetic systems. We are aware of the fact that our 
models do not simulate spin dynamics. However, any spin spin interaction as for example 



5 



in the incommensurate AFM structure of FeSb 2 S 4 has to be expected between the states 
given by the AFM, FM and NM models. Considering the orthorhombic structure which 
has four MnSb 2 S 4 formula units, two AFM configurations were accounted for, i.e. with 
the spin aligned oppositely in MnS@ octahedral chains, this will be called hereafter AFM1 
and another one with spins aligned parallel within a chain and oppositely between chains 
(AFM2). As for the monoclinic variety, the unit cell has two different manganese sites Mnl 
and Mn2. This leads to a first possibility which is to account for antiparallel spin alignment 
between Mnl and Mn2 sites leading to a ferrimagnetic (FIM) order. The other possibility 
is to double the cell along the third lattice vector c with Mnl and Mn2 all up-| in the first 
cell and Mnl and Mn2 all down-j in second cell, i.e. conforming with the spin spiral found 
for MnSb 2 S4 (mC28) 31 . Needless to say that the symmetry breaking due to the magnetic 
lattice orderings among Mn( j) and MnQ) in both structures computations are much heavier 
to carry out whence the limitation in the Brillouin zone integration in k points presented 
above. 

Further information about the nature of the interaction between atomic constituents can 
be provided using overlap population (OP) leading to the so-called COOP (crystal orbital 
overlap population) 32 or alternatively introducing the Hamiltonian based population COHP 
(crystal orbital Hamiltonian population) 33 . Both approaches lead to a qualitative description 
of the chemical interactions between two atomic species by assigning a bonding, non-bonding 
or antibonding character. A slight refinement of the COHP was recently proposed in form 
of the "covalent bond energy" E C ov which combines both COHP and COOP so as to make 
the resulting quantity independent of the choice of the zero of potential 34 . The Ecov was 
recently implemented within the ASW method 35 . Our experience with both COOP and 
Ecov shows that they give similar general trends although COOP exagerate the magnitude 
of antibonding states. We shall be using the E C ov description of the chemical bonding. 

IV. CALCULATION RESULTS AND DISCUSSION 
A. Total energy and magnetic moments 

Charge transfer is observed from Mn towards Sb, S and the empty spheres; nonetheless 
its amount is not significant in terms of an ionic description (such as Mn 2+ ), which is rarely 
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observed in the framework of such calculations. A more meaningful picture is provided from 
the quantum mixing of the valence states as it will be shown in the plots of the density 
of states (DOS) and the chemical bonding (E C ov) in nex t sections. The two polymorphs 
show similar trends concerning the total energy calculated for the non magnetic (NM) and 
spin polarized (SP) ferromagnetic -FM- and antiferromagnetic -AFM- models. Further fer- 
rimagnetic calculations in the monoclinic system were carried out. This is detailed in table 
1 which presents the results obtained after self consistent computations for the different 
magnetic configurations considered. For both modifications the FM state is favored with 
respect to the non magnetic one. The large gain in energy arises from the magnetic exchange 
of coupled high spin Mn 2+ when spin polarization is accounted for. In FM configuration 
the resulting total magnetization per formula unit is close to 5 For formally Mn 2+ two 
configurations are possible for the spin arrangements within the octahedral field of sulphur: 
A high spin -HS- configuration t% , e 2 with 5 unpaired spins and a low spin -LS- one: t\ , 
e° g resulting in only one unpaired spin. From this it can be suggested that at least in the 
ferromagnetic state divalent manganese is HS. However the total magnetization arises from 
all constituents of the lattice (Table 1) with the main contribution due to manganese. The 
SP-FM Mn moment amounts to 4.38 [i B within (oP28) MnSb 2 S 4 and 4.23 fi B for (mC28) 
MnSb2S4. This agrees with the value calculated for Mn 2+ in MnTe 21 and lies in the range 
of further LSDA results on high spin Mn 2+ from 4.3 to 4.6 {ib 21 ~ 23 - When the ferrimag- 
netic configuration is accounted for within the monoclinic variety a further stabilization is 
obtained and there is a cancelling out between moments so that total magnetization is zero. 
The AFM configurations show a further energy lowering for both varieties with a smaller en- 
ergy difference in the orthorhombic structure. This would suggest a lower Neel temperature 
for the high pressure orthorhombic variety. Within this structure the AFM1 configuration 
with the spins aligned oppositely in MnS 6 octahedral chains is found to be favored with 
respect to the AFM2 one (i.e. with spins aligned parallel within a chain and oppositely 
between chains). The small lowering in the moment carried by Mn 2+ (4.34 /j, b ) agrees with 
the value observed for a-MnS (4.1 Hb) 22 - For MnSb 2 S 4 (mC28) we find the same order of 
energies. Thus, the applied method predicts the preference of an antiparallel coupling of the 
spins of Mn-d electrons in a high spin state against a parallel coupling. The close magni- 
tudes of the moments between the FM and the AFM configurations lead to propose that the 
magnetic order might be Heisenberg-like whereby the magnetic susceptibility should obey a 
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Bonner-Fischer behavior 36 which is characteristic to linear spin chains. 

B. Non-spin polarized calculation DOS and chemical bonding 

The suggested NSP situation for MnSb 2 S 4 (oP28) and (mC28) results in a metallic 
behavior, analogous to studies on MnS and MnS2 20 ~ 22 . The site projected DOS are shown 
in fig. 2 a and b. The highest occupied states cross the Fermi level E F at a high density of 
states which is attributed to t 2g states from a crystal field analysis of Mn d states projections 
given in fig. 2c. These Mn 2+ t 2g states are only partly occupied by five electrons. The next 
bands above E F are formed by the Mn e 9 -states (fig. 2c). The splitting of the e g states results 
from deviations of the MnS6-groups from octahedral symmetry. Antimony and sulphur p- 
states form broad bonding states with the metal states in the energy range [-6,- leV] (see 
next paragraph). The DOS at low energies are s-bands of Sb (-10 eV) and S (-15 eV); the 
latter are found at lower energy due to the higher electronegativity of sulphur as compared 
to antimony. 

The chemical bonding within both orthorhombic and monoclinic MnSb2S4 are examined 
in the framework of the Ecw 34,35 for Mn-S, Mn-Sb and Sb-S pair interactions. The corre- 
sponding covalent bond energy Eqov plots are given in fig. 3 a and b. Negative, positive 
and nil Eqov magnitudes are relevant to bonding, antibonding and non-bonding character- 
istics. From this the major part of the valence band VB is bonding due to Mn-S interactions 
as well as to Sb-S albeit with a smaller magnitude; this contributes to the stabilization of 
the crystal lattice. The Sb-S interaction is observed with smaller magnitude and it remains 
bonding within the conduction band above Ep. This somehow provides an illustration for 
the description of the bonding given in the crystal structure section above. Mn-Sb interac- 
tion plays little role -as with respect to Mn-S one- within the major range of the VB. At 
the top of the VB the system becomes largely destabilized as the Fermi level is reached, 
i.e. where a large Mn-S E C ov as we h as Mn-Sb antibonding interactions with smaller mag- 
nitude can be observed. Although a large part of the Mn(t 2j ) are not engaged into Mn-S 
antibonding interaction in asfar as they are responsible for the onset of the Mn magnetic 
moment, the non magnetic configuration is clearly not favored from that. Lastly Mn-Mn 
interaction were observed too but with a much smaller magnitudes than all other explicited 
ones in both crystal varieties, so they are not shown here. Nevertheless it will be discussed 
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below that these bonds can have consequences on the electronic structure (cf. section III.C 
particularly for the monoclinic band structures). 

C. The electronic structure of spin polarized MnSb2S4 

1. Ferromagnetic state 

As shown by the site projected DOS in figs. 4 a and b, the spin polarization causes 
Mn 3d levels to split into majority spin (j) states which are lowered in energy relative to 
minority spin (J.) states at higher energy. Majority Mn d spin states completely lie below 
E F , thus being fully occupied by 5 electrons. The minority Mn d states are found above E F 
thus being completely empty. This indicates a closely non metallic situation with a small 
energy gap in the orthorhombic variety which reduces to a closing in (mC28) MnSb 2 S 4 . The 
DOS for manganese in both varieties exhibit peaks which closely resemble the ^(D-e^t) 
manifolds. Thus, the highest occupied states in the valence band are formed by Mn up 
spin e g states and the lowest unoccupied ones by down spin t 2g states. Concerning Sb and 
S DOS the latter can be observed to closely follow the shape of Mn pointing to the Mn-S 
coordination, i.e. with MnS 6 octahedra within which the major part of the bonding within 
the lattice occurs as discussed above. Spin polarization mainly affects Mn states so that 
there is hardly any energy shift between ( j) and (j) spin populations for Sb and S although 
residual moments were computed in both orthorhombic and monoclinic systems (Table 1). 

2. Ferrimagnetic (FIM) model in MnSb 2 S 4: (mC28) 

A first possibility to account for antiparallel spin alignment within (mC28) MnSb 2 S 4 was 
to allow for it between the two singly occupied Mn sublattices within the base centered 
monoclinic structure. The resulting energy differences shown in Table 1 are found in favor 
of this FIM configuration by 13.3 meV with respect to FM. The magnitudes of the moments 
are within range of FM calculations but the resulting magnetization is zero. The DOS and 
band structure given in fig. 5 a and b respectively. The DOS plot shows some similar 
features to FM (fig. 4b) but there is now a gap opening in the minority spins whereas 
a metallic behaviour is observed for majority spins. From the band structure plot in the 
same energy window the gap of ~0.6 eV can be observed between the VB and the CB in 
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U(A-E) direction which is along the k z axis of the monoclinic Brillouin zone. It is along this 
direction that the metallic behavior is obtained too as resulting from the crossing of single 
bands from the VB and the CB due to a large dispersion. Thus the monoclinic system, in 
an intermediate magnetic state (see relative energies in Table 1), is not a semiconductor but 
a half- metallic ferrimagnet with a relatively low DOS at E F due to single band crossing. 

3. Antiferromagnetic (AFM) models 

For all systems the energy differences shown in Table 1 are in favor of AFM ground state 
configurations (AFM1 for oP28). The result of enforced AF configuration is that the total 
up spin and down spin projected densities of states present the same contributions. As a 
consequence plots for one magnetic sublattice within each structure will be shown. In a 
narrow energy window around the Fermi level meant to exhibit the relevant features of the 
AFM ground state, figs. 6 and 7 give the DOS and band structure for orthorhombic and 
monoclinic AFM MnSb 2 S 4 respectively. The MnSb 2 S 4 (oP28) projected DOS (fig. 6) show 
a larger splitting around E F than in the FM DOS (fig. 4a). The larger gap is likely to arise 
from a shift of unoccupied minority Mn states to higher energies within the CB which can 
be a result of Mn-Mn interactions throughout the MnS@ chains. From fig. 6b showing the 
band structure its magnitude amounts to ~ 0.7 eV between T VB and Y C b for instance in 
the orthorhombic Brillouin zone. This results in a nonconducting state. Note that this gap 
for the AFM state is close to the experimental value of 0.77 eV 20 . Our calculations indicate 
the preference of an AFM configuration (AFM1, cf. energy differences in table 1) based on 
a simple model of alternating Mn moments along the rods. This is somehow similar to the 
a-MnS case examined by Tappero et al. 23 . 

AFM ground state site projected DOS of MnSb 2 S 4 (mC£5) (fig. 7) show different features 
from the ferrimagnetic case (fig. 5a) because both Mnl and Mn2 are now polarized up or 
down within a magnetic sublattice (see for instance the change of orientation of Mnl and 
Mn2 DOS above E F ); this results in larger n(E F ). In termes of band structure (fig. 7b) this 
involves enhanced band crossing along AE direction (along k z direction) as it can be observed 
from the confrontation with the ferrimagnetic band structure (fig. 5b). From such a band 
dispersion and crossing the system is obtained as weakly metallic. This is somehow opposed 
to the semiconducting state proposed experimentally. Nevertheless both monoclinic and 
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orthorhombic varieties have been shown to possess similar features and the final answer on 
the question for the coupling of the magnetic moments will be given by neutron diffraction. 
Related investigations are in progress 31 . 



V. CONCLUSION 



The electronic structure of MnSb2S4 in both, the orthorhombic and the monoclinic mod- 
ifications were calculated within the local spin approximation for non magnetic as well as 
for spinpolarized ferromagnetic, ferrimagnetic and antiferromagnetic models. According to 
total energy calculations the spin polarized states with high spin Mn 2+ are largely preferred 
to a non spin polarized one (Table 1). Magnetic moments of ~4.3 hb are calculated in agree- 
ment with high spin Mn 2+ configuration known from MnS and MnS2- For both MnSb 2 S4 
varieties the AFM model shows an additional energy gain, thus becoming the ground state. 
These results are accompanied by significant differences in the electronic structures of the 
models. The NSP model leads to a metallic behavior for both modifications with a partly 
filled VB formed by Mn t 2g and the CB by the empty Mn e g bands shown by a crystal 
field analysis. In the orthorhombic system FM and AFM models lead to the experimentally 
observed semiconducting characteristics with a larger gap obtained for the AFM ground 
state. Differences in the electronic structures concerning the CB and the VB are due to the 
crystal structures. Calculations for MnSb 2 S4 (oP28) reveal a band gap of 0.7 eV, close to 
the experimental value of 0.77 eV. In MnSb 2 S 4 (mC2#) two Mn sites are present which have 
a significantly different environment by sulphur and therefore the site projected DOS for 
Mn shows a broadening, hence the VB is broadened too in comparison to the orthorhombic 
modification. On the other hand, the empty minority spin Mn d states in the conduction 
band are sharper for the monoclinic modification. This is related to the higher local symme- 
try at the Mn sites. The computed intermediate ferrimagnetic state exhibits a half metallic 
behavior due to single Mn bands crossing along the AE direction in the Brillouin zone, i.e. 
along k z . This is enhanced in the AFM ground state. Although the antiferromagnetic na- 
ture of the ground state of both modifications of MnSb 2 S 4 becomes evident by the present 
calculations, further investigations of electrical conductivity to reveal the semiconducting 
properties are needed, they are underway. 
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TABLE I: Crystal data from literature and calculation results for orthorhombic (oP28) and 
monoclinic (mC28) MnSb2S4 -NSP=non spin polarized; FM=ferromagnetic; FIM=ferrimagnetic; 
AFM= antiferromagnetic. 
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FIG. 2: a) and b) Show the site projected DOS for one formula unit of non magnetic MnSb2S4 
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FIG. 3: Chemical bonding properties from covalent bond energy J^cov approach within MnSb2S4 
per formula unit: a) oP28 orthorhombic variety, b) mC28 monoclinic variety for one of the two 
manganese sites -Sb and S regroup partial contributions from all lattice sites. 
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FIG. 4: Spin resolved site projected DOS per formula unit for a) MnSb2S4 (oP28), b) MnSb2S4 
(mC28) -Sb and S regroup partial contributions from all lattice sites. 
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FIG. 5: a) DOS (Ml and M2 stand for Mnl and Mn2 respectively; Sb and S regroup partial 
contributions from all lattice sites) and b) band structure in a narrow energy window around the 
Fermi level of ferrimagnetic intermediate state of monoclinic MnSb2S^(mC28) (Solid lines (|), 
dotted lines (].)). 
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(b) 

FIG. 6: DOS and band structure in a narrow energy window around of antiferromagnetic 
ground state of orthorhombic MnSb2S4 (oP28). 
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FIG. 7: DOS and band structure in a narrow energy window around of antiferromagnetic 
ground state of monoclinic MnSb2S4(mC2#). 
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